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Porque não dar uma vista de olhos ao nosso pla- 
neta a partir de outro ângulo? 

Para tal é preciso ir onde haja explorações diárias 
dos polos ao fundo do mar. E isso só no Edificio 
da ITT em Nova lorque. 

Aí a ITT prossegue nessa tarefa com a colaboração 
das suas fábricas espalhadas em mais de 50 países 
do mundo, nas quais mais de 150 000 pessoas dão 
o melhor do seu esforço para conseguir novas 
técnicas e novos processos. Mas que género de 
trabalho é esse? São as telecomunicações e a elec- 
trónica. Na ânsia de um mundo melhor e unido 
os homens procuram ligar-se, entender-se, comuni- 
car. E as companhias da ITT fornecem esse meio. 
Conquanto que cada fábrica produza vários tipos 
de materiais, há, todavia, um ou outro produto 
especial que abastece todas as fábricas do Sis- 
tema. 

Assim, em Lisboa fabricam-se radiossondas e emis- 


sores de comunicações SSB; em Madrid, cabos e 
telefones; em Paris, equipamentos de telecomuni- 
cações troposféricas telefónicas; em 
Viena, intercomunicadores; em Antuérpia, sistemas 
de tubos pneumáticos para rápido transporte de 
correspondência; em Copenhague, equipamentos 
radiomarítimos; em Estugarda, Alemanha, equipa- 
mentos para navegação aérea e transmissores; em 
Milão, equipamentos de microssondas; em Oslo, 
frigoríficos, fios e cabos; em Estocolmo, equipa- 
mentos telefónicos variados; em Zurique, compo- 
nentes: em Croydon, Inglaterra, teleimpressoras; 
em Londres, cabos submarinos e receptores de 
rádio e TY; em Chicago, equipamento militar; 
em Buenos Aires, rádio e TY; no Rio de Janeiro, 
electro-domésticos; na cidade do México e no 
Chile, aparelhos telefónicos, etc. Muitas outras 
companhias estão espalhadas pela África, Japão, 
Austrália, etc. 
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Peças em borracha podem agora ser produzidas em equi- 
pamento convencional de processar plásticos com resinas 
copoliímeras ALATHON E/VA de etileno e acetato de vinilo. 
As resinas ALATHON E/VA possuem a resiliência e a 
resistência à flexão repetida da borracha e de outros ma- 
teriais plásticos flexíveis sem quaisquer problemas de migração, 
pois não contém nenhum plastificante. 


As resinas ALATHON E/VA são também fáceis de mol- 


novos mercados lucrativos 


dar por sopro, de injectar, de exturdir ou de moldar por 
vácuo. À rapidez dos ciclos nestas operações oferece ainda 
duas vantagens importantes : economia de tempo e custo. 
Outras características : tenacidade a baixas temperaturas, 
fácil aceitação de cargas e fácil mistura com outros mate- 
riais termoplásticos. 

O ALATHON E/VA tem um lugar importante na vossa 
gama de materiais de moldação. 
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TERMOESTÁTICA 
por A. GOUVEA PORTELA 
Prol, L:S. T. 
RESUMO SYNOPSIS 
Este artigo reformula a quiomálica da Termoestática e This article reformulates the axiomalique of the Ther- 
introduz a aplicação dos métodos de Programação (nomea- mostatic and introdures the application of Programation 
damente linear) a resolução de problemas termoestáticos. Methods (ex. grat, linear Programation) to the solution of 


Thermostatics Problems. 
0 — INTRODUÇÃO 


Como toda a Teoria Física, a Termoestática consente várias axiomátiças. 

Nem sempre a axiomática proposta inicialmente é a mais conveniente e, com o decorrer do 
tempo, surgem outras formas de apresentar as bases duma teoria que conferem a esta um renovado 
interesse. 

O objectivo desta nota é justamente sugerir uma axiomática para a Termoestática que parece 
possuir certas vantagens que submetemos ao julgamento dos interessados nestas questões. 

Por maior consistência que se deseje dar a uma teoria, há sempre um certo número de termos 
cujo sentido subjectivo se supõe conhecido, assim vamos supor dominados os conceitos correspon- 
dentes a vocábulos, tais como: universo, espaço no sentido Topológico, contorno ou fronteira, 
conexão, referencial, estacionário, etc. 


1 APRESENTAÇÃO DA TEORIA 


Sistema Termoestático é toda a parcela finita do Universo delimitada por um contorno fechado 
constituindo este uma figura geométrica conexa no espaço Euclideano e satisfazendo às seguintes 
propriedades : 


a) A fronteira tem propriedades invariantes no tempo. 

b) O Sistema está contido nessa fronteira há co tempo. 

c) O Sistema (5) é por seu turno uma topologia (anel) com as seguintes propriedades 
estruturais internas (organização): 


5 é mensurável e tem pelo menos (n +-1) medidas finitas =, todas completamente aditivas. 
d) Existe um número n inteiro e finito de medidas independentes p; (i=1,...n) que dum 
modo completo caracterizam o sistema S. 
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e) Entre as medidas » de S existe uma também completamente aditiva que designaremos por 
energia interna U que tem a propriedade de tomar o valor mínimo compatível com as condições 
impostas pela fronteira. U é um escalar. 


Embora a axiomática de um sistema lógico não necessite formalmente de justificação, contudo, 
para que tenha interesse físico, é usual e conveniente apresentar uma interpretação ou tradução em 
metalinguagem a fim de mostrar que há na Natureza sistemas reais cujo comportamento se apro- 
xima razoávelmente do modelo formal apresentado. 

Neste sentido apresentamos, a seguir, esclarecimentos e justificações dos axiomas escolhidos : 


— As proposições a) e b) garantem que o sistema Termoestático isolado na fronteira atingiu um 
estado estável no sentido macroscópico. 

Deste modo todos os fenómenos de relaxação e de histeresis se supõem completados e a inva- 
riância das propriedades da fronteira elimina o ruído do universo circundante que, para todos os 
efeitos, tem as propriedades da fronteira. 


— À proposição c) é fundamental porque postula a existência de medidas para S dum certo tipo 
e caracteriza essas medidas com a propriedade da sua somabilidade. Isto é: 


3 
Se A e B forem dois conjuntos satisfazendo a: do R 
AnB=O 
CLAUB:S 


então v(AUB)=vê(A) + (B) 
e ainda 
Se for dado uma série numerável de conjuntos, Er, disjuntos dois a dois então: 


p 2 E=2Zu (E) 
r Fr 


com E :. S e PAP 
Er n e = 6) r == k 
E, U Ex U Es is E 5 


— À Proposição d) tem o mérito de tornar finito o número de medidas macroscópicas a efectuar 
no sistema termoestático para o caracterizar completamente sob o ponto de vista termoestático. 

Note-se que se não indica quantas são essas medidas. 

O número de medidas dependerá do problema a resolver e da precisão desejada e caberá ao 
físico escolher essas medidas; unicamente se postula que esse número é finito. 


As proposições referidas até aqui (a, b,c e d) permitem extrair algumas consequências que 
que convém explorar desde já para boa inteligência do que se segue. 

Qualquer medida go de S satisfazendo a c) que não esteja incluída no número das n medidas 
v; consideradas independentes é uma função das n medidas gi 


o = Po (ma... Un) 


Tendo em vista a proposição d), vo será uma função homogénea de 1.2 ordem de 4... tn,e 
continua no ponto considerado se o for completamente aditivo. 

Esta propriedade resulta de se ter postulado que as medidas » são completamente aditivas, o 
que implica que o conjunto dos pontos de descontinuidade tenha medida nula. 
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A partir das medidas vo, ti ...vn todas completamente aditivas, sendo n independentes, 
pode construir-se um espaço de representação En,1 (eixos rectilíneos ortogonais, centrado e afim). 

Nesse espaço pode ser representada a hipersuperfície vo =P (H1...tn) onde estarão 
situados todos os pontos que correspondem aos estados atingíveis pelo sistema termoestático. 

Também nesse referêncial será possível representar quaisquer constrangimentos de forma Ok (po, 
VI, cc En) E O. 

Estes hiperespaços de constrangimento da forma 9x (po, t4,+-- Hu) = O interceptarão a hiper- 
sufície Lo == ?o (14 ... Hu), restringindo o domínio de existência dos pontos representativos do estado 
do sistema. 

A proposição e) tem a função de fornecer um critério para verificar se determinado estado 
representado por um ponto no referencial En é ou não estável. 

Assim o ponto P satisfazendo à hipersufície U =U (p4... tu) e a todos os contrangimentos 
de forma Ok (2,5, Hi... Yu) = O, e que se admitem compatíveis para que P tenha existência, não fica 
bem determinado, em geral. 

Com efeito, a menos que os constrangimentos definam o ponto duma forma unívoca, em geral, 
o que se verifica é a existência de um domínio e só um critério adicional permitirá escolher um ponto 
nesse domínio. 

A função da proposição e) é precisamente fornecer esse critério, o qual se formula desta forma 
simples: U será mínimo. Se só existir um ponto que minimize U e for esse o ponto que representa 
o estado do sistema então esse estado é estável. 

Em resumo, toda a termoestática se reduz a um problema de max e min condicionados, ou 
seja a um problema de programação não linear, o qual consiste em procurar num domínio delimitado 
um ponto que optimise determinada função desse ponto. 


2 — EXPLORAÇÃO DO SISTEMA LÓGICO 


Uma vez que a axiomática da teoria foi apresentado e sumáriamente interpretada, convém pas- 
sar à sua exploração. 
Para o efeito abordaremos três temas 


Conceito de matriz de rigidez 
Espaço 2n 41 
Conceito de entropia 


2.1 — Conceito de matriz de rigidez 


Embora o referencial de representação seja um E, (ortogonal, centrado), nem por isso deixa 
de ter interesse caracterizar a variância das grandezas intervenientes. 

Assim todos os »; são ou escalares ou afinores contravariantes. 

A energia interna U é um escalar e tem as dimensões de uma energia. Porque U é uma fun- 
ção contínua de |'4... !y será ainda: 


- d U A . — Á A / 

designaremos por 7; = a =; será ainda uma função contínua de 4... tn, porém 7; é uma função 
Hi 

homogénea de grau zero pois é uma relação de duas medidas que gozam da propriedade de somabi- 
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lidade, porém =; não goza dessa propriedade. Se designarmos as variáveis | de extensivas, as variá- 
veis x; tomarão o nome de intensivas. A expressão de dU pode ainda escrever-se 


dU = [ri]. [dpi] 
(1, nm) (n,1) 
Porque 7; é uma função de pi ...pn será 
Toi " Jé 
d tj == 2 Edi: . d uk == b> é. MO . d px 
E O pk k dpi dp 
e 
Dim 92 U - 
[d =] =| SE | tam =| 55 a Ci 

O pa d pi O pk 


(n, 1) (n,n) (n,1) 


d E 4 . 
A será ainda uma função continua de mM... Pu. | 
e Em 


d Ti 


| tem o nome de matriz de rigidez. 
Ô pk 


“0x 
Esta matriz rigidez | 


d dk 
das variáveis dz; comas dri. 
Em cada ponto, representativo dum estado estável, a matriz de rigidez converte-se numa 


matriz numérica. 


| pode ser invertida parcial ou totalmente, permutando a posição 


* a 6 
O Hj O pk dvk dO pi 


A circunstância de U ser contínuo arrasta a igualdade : e daí re- 


sulta que a matriz de rigidez é simétrica. 
“Sm Tt 
Finalmente. | tem a designação de matriz de elasticidade. 
Pk 
Note-se que toda a teoria de elasticidade se propõe por meio de uma matriz de rigidez 
de forma 
[T] = [MJ]. [5] 
onde 
[5] representa o tensor de deformação 


[T] representa o vector de deformação 
[M] representa a matriz de rigidez 


2.2 — Espaço Ex. 


A introdução das variáveis 7; (intensivas) vai permitir a representação do sistema termoestá- 
tico S num referencial mais amplo com 2n + 1 dimensões. 

As coordenadas serão agora U, gi ... tn, 1... Tn. O hiperespaço termoestático correspondente 
à Superfície Termoestática será representado neste referencial por: 


U=U (moto, Ta) nt 1 
e TE== Ty (Ugo Em; Ri seo Fa ) equaçoes 
E 


Os constrangimentos, em número de m, tomam a forma seguinte; 


dr lisa Ei, Tita) E O 
K=1, em 
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Note-se que o número de graus de liberdade do Sistema Termoestático não constrangido con- 
tinua a ser n= (2n + 1) — (n + 1), como não podia deixar de se verificar uma vez que se postu- 
lou ser n o número de variáveis independentes. 

De novo o problema da estabilidade dos sistemas termoestáticos se formula extremando 
(minimizando) a função U, condicionado a n constrangimentos =; e m constrangimentos Sk, o que 
corresponde a um problema de max. e min. condicionados ou ainda a uma programação não linear. 


2.3 — Conceito de entropia 


Em termodinâmica (irreversível) a entropia desempenha uma função primordial. 

Contudo na axiomática da termoestática o papel da entropia é mais modesto. 

Com efeito, a circunstância de se ter imposto que o número de medidas a efectuar ao sistema 
termodinâmico deve ser finito (n) resulta que é necessário apurar quais as medidas |; (variáveis 
extensíveis) que influem no sistema o que pode experimentalmente determinar-se procurando esta- 


belecer as diferenciais RE um Fi. 
Hi 

A natureza do problema vai permitir arrolar assim quais as variáveis =; que possuem uma 
correlação significativa com U. 

Mas a circunstância de ser práticamente co (e não finito como se postulou) o número de 
variáveis influentes em U é formalmente necessário acrescentar mais uma medida (que designare- 
mos por |n) que representa a totalidade da interacção restante, não expressamente contida nos 
parâmetros pi... Un-l. AW 

Assim haverá (n—1) medidas »; môgíoscópicas mensuráveis experimentalmente e mais uma 
medida vn não mensurável directamente à qual corresponderá uma propriedade intensiva 
dU 


d Em 


e que é responsável pelo trabalho rn dyn que representa a energia não contabilizada no 


somatório 2 mi; dpi. 


A esta medida un dá-se o nome de entropia, a Tn O de temperatura e a Tn dtn, com a grandeza 
de uma energia, a designação de calor. 

À circunstância de ser possível conceber deslocamentos do Ponto P (pi... vn) a dyn=0, 
aos quais correspondem energias rn dyn = 0, não invalida a asserção feita, porque para o sistema 
lógico ser consistente é necessário que seja aplicável mesmo que din -fO. 

Dada a liberdade de escolha que dispomos quanto a yy e tn poderemos impor que 7; > 0. 

Em resumo, poder-se-á escrever: 


n—1 ) o 
dU = > eu d ui + Lá dun 
j=1 O pi O Um 
ou numa simbologia mais corrente 
n—1) 
dU= 3 edi + TAS 
i=10pi 


Este modo de proceder é frequente em física. Um número suficiente de experiências foram feitas 
sobre um sistema dado, foram oferecidas, à priori, um certo número de variáveis consideradas 
influentes no sistema, procuraram-se correlações entre a função (U) e essas variáveis 1»; e assim se 


* * . A * d U . 
determinaram os coeficientes de influência —— , (i=1,...n—1). 
dpi 
Porém ficou no fim um resto. Então é legítimo oferecer mais uma variável |, que já se não 


pode medir directamente mas que é tornada responsável pela interacção restante. 
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Deste modo, é possível representar o estado do sistema num referencial com (n + 1) dimen- 
sões, por um ponto situado numa superfície bem determinada. 

Para concluir tem interesse chamar a atenção de que, nos sistemas termoestáticos no estado 
estável, a circunstância de U ser mínimo implica ser »,=S máximo. 

Com efeito, 


nl 
> id representa a energia correspondente à totalidade dos parâmetros extensivos men- 
j= f ” ” F ” 2. 
suráveis experimentalmente e é bem determinada para cada deslocamento virtual 
)P do ponto P (»...n) situado na superfície termoestática e sujeito aos con- 
trangimentos impostos. 
n—l 
Donde CU—T/S = 2 t;9g; e o segundo membro é bem determinado para cada desloca- 
| 


mento )P e porque T >0O (por hipótese), e)U <'0 (seu for minimo), terá de ser 95 < 0 (ou S 
máximo), conforme se considerou U = U (g... vn-1,S)ouS=5S (ui... py U). 


3 — APLICAÇÃO DE PROGRAMAÇÃO LINEAR 


Nos capítulos anteriores foi mostrado como os problemas de termoestática se reduzem ao da 
programação não linear. Infelizmente não há um método universal de resolver problemas deste tipo 
mas em contrapartida há métodos perfeitamente explorados para tratar a «programação linear». 
Vejamos em que condições é possível e legítimo aplicar este segundo método. 

Se a vizinhança explorada em torno do ponto (no espaço de representação) é relativamente 
reduzida é, em geral, legitimo: 


— linearizar U=o0 (ij. En, Tio Tn) 
E na é RE 
— considerar a matriz ei como tendo elementos constantes. 
VXj 
' " à | " > 
— linearizar o sistema de equações: Pe (141 cestins RL cos Tn) ( Ss) 0 


Então, a aplicação da programação linear é possível e não oferece dificuldades e propõe-se da 
forma seguinte : 


Sendo dado o sistema de constrangimentos lineares seguinte: 


o = | |. [op + | SE Lóm] 


O Fi |] Ri. 


minimizar a função critério (também linear) seguinte : 


“q : 7 
U=[ | Dul+| So). tal 


| du d Ti 


LP] 


Sendo todos os deslocamentos a partir de um Ponto P de equilíbrio tomado para ponto 
de partida. 

Note-se que antes de aplicar o método ter-se-ia que averiguar se o poliedro formado por todos 
os hiperplanos era convexo, bem como outros cuidados formais, que aqui não nos referimos 
expressamente. 
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4 — INTERESSE DO MÉTODO 


As condições de fronteira que era usual impor aos sistemas termoestáticos consistiam tradi- 
cionalmente nos seguintes tipos : 


ou dt = 
ou dm =o0 


Isto é, ou se impunham a P deslocamentos a volume, entropia, componente (i), constantes 
ou a pressão, temperatura, potencial químico, constantes. 


Nestas condições não há interesse em recorrer ao método de programação linear, pois que 
bastava estudar a matriz de rigidez ou a sua inversa (parcial ou total) para concluir da estabili- 
dade do sistema. 

É também clássico o emprego da transformada de Lagrange-Legendre da função U ou S, tais 
como a entalpia, função Gibbs, funções de Massieu, para estudar estes casos, que correspondem a 
inverter parcial ou totalmente as matrizes de rigidez ou de elasticidade. 

Mas se as condições de fronteira se apresentarem sob formas mais complexas, como: 


SU-Zp=5õ. 
a t+ KkRP-—T< O 
P>3 


Então só o formalismo da programação linear tem potência para resolver esta classe de pro- 
blemas. 

Pergunta-se, mas haverá casos reais em que tal suceda? 

A automação com as suas múltiplas e variadas aplicações ao controle de processos e máqui- 
nas impõe condições que já se não descrevem simplesmente por formas tais como d p= 0 ou 
d =; = 0, então tem completo e cabal cabimento o emprego de programação linear. 


REFERÊNCIAS: Measure Theory — Halmos. 
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C. D. U. 624.314.5/54 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL ') 


I — Breve nota mensal 


Durante as 3 semanas iniciais de Setembro prática- 
mente não choveu, pelo que as afluências desse pe- 
ríodo vieram a reflectir, como seria de esperar, o natu- 
ral prolongamento do regime estival anterior, o qual 
já apresentava carácter de extrema secura. Na última 
semana de Setembro, porém, a ocorrência de apreciá- 
veis precipitações ocasionou uma subida de afluências 
de certo modo sensível, sobretudo nos últimos dias 
do mês, 

As afluéências dos primeiros 20 dias de Setembro 
(apenas 48 GWh) correspondem a um regime em tudo 
idêntico ao de Agosto (73 GWh). Nos últimos 10 dias 
do mês as afluências atingiram 81 GWh dos quais cerca 
de so GWh concentrados nos últimos 3 dias. A afluén- 
cia global do mês de Setembro perdeu assim o seu 
carácter extremo muito embora se tivesse mantido 
ainda na zona dos regimes desfavoráveis. 

O coeficiente de produtibilidade hidroeléctrica foi 
apenas de o,7o, tendo as respectivas afluências uma 
probabilidade de serem excedidas da ordem dos 91º/,. 


II — Elementos gerais 
Acumulados desde 1 de Janeiro 


GWh Variação 
1964 | 1965 | 9% 

Produção hidráulica (Ph)... 32445 25610 — 21 
Produção térmica (P+).....| 111,8) 380,9 | + 242 
Produção total (PT). ..... 3955,8| 29424 — 12 
Energia recebida de empresas " , 
não pertencentes ao RNC (Er) | 20,9 15,6 — 24 
Exportações (Ex)... ...... 36,5 34, — 91 
Importações (1) ......... 0,1 439,8 | — 
Saldo importador (St). ..... 136,4 4+435,9 — 
Consumo em bombagem (Cp) ') re 6, | ESA 


Produção para con- o : À 
sumos perman. (Pep) (2) 2679,4| 3046,8 | + 15,1 


Produção para con- 


sumos não perman. (Pc 3) | 660,5, 340,3] -— 48 
Tot caes nine ii 3399,9 | 3993,9 + 2 
Coeficiente de hidrauticidade | 112] 0,68 — 


NOTAS 

(1) Energia consumida em bombagem para armazenamento na Al- 
bufeira do Alto Rabagão. 

(2) A produção para consumos permanentes (Pcp) é determinada 


Pep = PT + Er + Si — Cb — Penp 


pela seguinte expressão : 


(3) Consideram-se consumos não permanentes, os indicados em 4.2. 


(4) O aumento percentual de produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, é 
respectivamente de 11,0ºJo e 14,3 o» 


II —Diagramas de carga dos dias característicos 


4 teira: 
1961 1965 
Produção hidráulica (Ph) MWh/ 11103 8954 
Produção térmica (br) MWh| 1030 2260 
Produção total (PT) MWh| 12133 11212 
Trocas com [ Export. (Ex) MWh U 15 
- Espanha | Import. (1) MWh 0 No 
Consum em bomb. hidroel, (C4) MWh Ú o qi 
Prod. para cons. perm.  (Lc,) MWh| 10847 12135 
Prod, para cons. não perm, (Penp) MWh | 1285 172 
TOTAL Pr (1x) MWh| 12153 12907 
do pra max, MW) 661 699 
Ee j *otência min,  MW| 319 333 
SE Pro Uritizda ponta horas | 18,4 18,5 
TE Factor de carga 0,76 0,77 
s — | Potênciamáx. MW| 608 | 670 
a p Potência min. MW 258 301 
Us e Utiliz. da ponta horas 18,0 18,1 
Factor de carga 0,75 0,75 


SETEMBRO 


2.3 “DIAGRAMAS DE CARGA DOS DIAS CARACTERÍSTICOS 


+ | 
7 dy ' EO a Emo 
: | R 
É) 


Pr elI-Ex)-C, 


HE 


1a 1Ix-6S5 
peléi M 
vovego VISÕES DS UT suas dO Ix-64 


e Produção porg consumos não permanentes (1965) 


IV —Energia armazenada nas principais albufeiras 


No fim do mês 
Albufeiras ; cdi SA 
GWh o (8) 
Alto Rabagão .. cc... “| 105,5 10,8 
Fama sis pads s%% z 0,0 0,0 
Venda Nova . c. cc. «| 928 121 
Salamonde . . .. «cc c«| 269 96,7 
Caniçádo co ccmo o o 31,1 94,0 
Vil susso» mn amo TRE 12,5 
RGE, é a é sa as é qo) MOS 32,6 
Castelo do Bode. +... .. 76,2 46,7 
CMMBONS é» é é mis é é 3,5 63,9 
Lagoa Comprida +. . +... 1,5 (2) 0,0 
Bimta LIA o cc us co. 34,0 55,2 
RN. & Dus ma sim as 3,4 26,4 
POVO sn uiicrãs 2,4 (3) 20,2 
Total com À, Rabagão ++ »| 500,5 23,1 
| sem À, Rabagão . . .| 3950 34,8 
NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras, 

(2) Inclui 1,3 GWh armazenados em Vale do Rossim no início do 
mês e 1,5 GWh no fim do mês, 

(3) Inclui 0,3 GWh armazenados no açude do Poio no início do mês 
e 0,3 GWh no fim do més. 


(*) Elementos extraidos das estatisticas mensais do Repar- 
tidor Nacional de Cargas (BR. N. C,), As produções e 
os consumos das empresas do R, N. O. representani 
cerca de 94 */, dos totais do Pais, 


Ponte rolante eléctrica de 100+-10 +t/14,2m, 
em serviço na central subterrânea de Alto Rabagão . 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 


TURBINAS HIDRÁULICAS —— Fabrico segundo licença de A. C.M. de Vevey, S.4. 
TURBINAS A VAPOR ——— —-— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR ————————— Projecto e fabrico segundo licença de Foster 


Wheeler, Co. 
EQUIPAMENTOS E INSTALAÇÕES 


INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGUE :... 


ALVERCÇCA DO RIBATEJO- PORTUGAL 


TECNICA XIH 


aços vazados 

ferros meehanite 
mecânica MÉDIA E PESADA 

construções metálicas 


COMPORTAS 
TURBINAS 


CONDUTAS FORÇADAS | 
RESERVATÓRIOS SOB PRESSÃO 
PERMUTADORES DE CALOR 


Indústria dos Petróleos. 


Coluna de 1000 mm 4 =< 10 000 mm, 
em fase final de fabrico. 
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COMPANHIA UNIÃO FABRIL 
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NOTA SOBRE OS DOMÍNIOS DE UTILIZAÇÃO DAS MÉDIAS 
TENSÕES NO TRANSPORTE DA ENERGIA ELÉCTRICA 


RESUMO 


OU artigo estabelece os doménios de utilização das médias 
tensões em função da potência a transportar e da distância 
de transporte. O diagrama apresentado tem interesse para 
uma definição prévia das tensões que podem interessar a 
um determinado transporte. 


1. Quando se procura definir qual a tensão 
adequada para um determinado transporte de 
energia é vantajoso efectuar um cálculo prévio 
que permita restringir o domínio que se investi- 
gará com pormenor. 

Para as médias tensões — por médias tensões 
entendem-se as que vão de 6 kV a 30 kV (!) — 
um dos critérios que se usa, para definir o domi- 
nio das tensões que interessam é o da queda de 
tensão: uma tensão será conveniente para o 
transporte se conduzir a quedas de tensão aceitá- 
veis (da ordem de 5 a 10“/o nas médias tensões) 
e isto sem obrigar ao uso de condutores com 
secção excessivamente grande. Adiante (alínea 3) 
se exporá o que se entende por secção excessi- 
mente grande. 

Para aplicar este critério, com comodidade, 
será vantajoso dispor de um diagrama que rela- 
cione as quedas de tensão, com a tensão de 
transporte Uc, a distância | e a potência activa 
transportada P. 

É certo que a escolha final da tensão terá de 
atender, entre outros factores, à temperatura dos 
condutores, ao custo da energia perdida na linha, 
e aos encargos com o capital a investir. Mas 


()— O artigo 83º do Decreto-Lei n.º 43335 impõe 
que as tensões a adoptar no transporte ou na grande dis- 
tribuição sejam 6000 V, 15 000 V, 30000 YV, 60000 V, 
100 000 V, 150 000 V, ou 220 000 V. 
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SYNOPSIS 


This paper establishes the field ofapplication of medium 
roltages as a function of the power to be transmitted and 
the distance of transmission, 

The diagram presented is useful in a previous definition 
of the voltages which can be involved in a power trans- 
mission. 


uma análise das quedas de tensão dá informação 
valiosa para uma apreciação prévia das condi- 
ções do transporte. 

Em comunicação apresentada nas «Jornadas de 
Engenharia de Moçambique» em Abril de 1965, 
pelo Prof. Eng.º António Carrisso e pelo autor 
deste artigo (?), relacionava-se o produto IP com 
a tensão de transporte Uc, usando a expressão 
(8), que adiante se deduz, e representava-se a 
função IP (Uc) para condutores com diferentes 
secções. 

Julga-se ter, também, interesse a representação 
da função P(l) para cada tensão, linearizando o 
traçado pelo uso de escalas logarítmicas. É esse 
o objecto desta nota, que se refere exclusiva- 
mente a linhas aéreas. 


2. Num transporte em média tensão, é, em 
geral, aceitável tratar a linha como um circuito 
com constantes concentradas e supor nula a sua 
admitância transversal. É o que se fará na aná- 
lise que segue. Igualmente se nota que se admi- 
tem sinusoidais as grandezas alternadas que 
intervêm no cálculo e apenas se consideram os 
regimes permanentes. 


(*) — António Carrisso e Domingos Moura — «Bases 
para o projecto de uma rede de distribuição rural», 
Comunicação n.º 40, Primeiras Jornadas de Engenharia 
de Moçambique, Lourenço Marques, 1965. 
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Pretende-se transportar uma potência electro- 
magnética constante P (watt) a uma distância | 
(metro) usando uma linha aérea trifásica equili- 
brada. Detinindo como «queda de tensão» U“ a 
diferença entre os valores eficazes da tensão sim- 
ples no princípio e no fim da linha, será, muito 
aproximadamente, como se sabe, 


Uo=rt1]Icosg+xlIsende... (1) 


em que 


U'o— valor eficaz da queda de tensão (volt). 

r — resistência unitária do condutor (ohm. 
metro” ?), 

— distância (metro). 


e 


| — valor eficaz da intensidade da corrente 
numa fase (ampere). 

é — argumento da impedância do receptor 
(radiano). 

x —reactância cíclica unitária da linha (ohm. 
metro !). 


Se for rc a resistividade do material condu- 
tor (ohm. metro) e s a secção do condutor duma 
fase (metro”), será 

Te 


f = saidas E dA 
: 


Notando que 


Pas Mol Dn nes s » « EB 
e que 
n 
mo cu umas ua é rica Rá 
100 
em que 
Us — valor eficaz da tensão simples, que ao 
aplicar as expressões (3) e (4) é suposta 
constante ao longo da linha (volt). 
m — queda de tensão (percentagem de Us). 


pode-se escrever 


m t IP KuLr 
Us emese 1 qe E tg À “ . . 5 
100 3 Us , 3 Us i ta) 


ou ainda 
m U2. 


IP o VER (é . “ ” (6 
100 (r + x tg 4) 
em que Uç é o valor eficaz da tensão composta 


(volt). 
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Usando as unidades 


| — quilómetro. 

P — quilowatt. 

r — ohm. quilómetro” !, 
x — ohm. quilómetro” |. 


Uc — quilovolt. 
a expressão (6) transforma-se em 


— I0m UZ% 
r-+x.tgó 


:P 


wo o 4) 


Nota-se que a expressão (7) supõe um receptor 
de potência P colocado no extremo da linha. Se 
a potência P estiver uniformemente distribuída 
ao longo da linha de comprimento 1 em vez da 
expressão (7) obtém-se a expressão 


20 m U* 


IP=———— 
r-+x tgó 


. » (8) 


Tomando como constantes a queda de tensão 
m, a reaçtância cíclica unitária x e o factor de 
potência da carga (cos &), pode escrever-se 


Es Ma mas Ega (9) 


À expressão (9) correspondem segmentos de 
recta no plano (], P), desde que as escalas de 1 
e de P, sejam logarítmicas. São esses segmentos 
de recta que se traçam na Figura 1, para a hipó- 
tese de cargas concentradas no fim da linha. 
Cada segmento corresponde a um condutor de 
resistência unitária r e a uma tensão Uc nor- 
malizada. 

O diagrama que se representa na Figura 1 foi 
traçado para condutores em cobre duro e su- 
pondo 


m = 5 
x = 0,4 ohm. quilómetro ! 
cos é = 0,8. 


Estes valores são os que habitualmente se 
consideram para uma análise sumária de um 
transporte em média tensão e correspondem às 
situações que, como é bem sabido, se apresen- 
tam com mais frequência e em que a expressão 
(1) conduz a uma aproximação aceitável. 


3. As faixas que na Figura 1 representam o 
domínio de cada tensão são limitadas por dois 
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segmentos: o inferior corresponde à menor sec- 
ção que se emprega em linhas de alta tensão; 
o superior procurou-se que coincidisse quanto 
possível com o segmento que, para a tensão 
imediatamente superior, corresponde à menor 
secção. 

A regulamentação portuguesa (!) proibe que 
em linhas aéreas de alta tensão se usem condu- 
tores em cobre com secção inferior a 10 mm?, 
Admite-se portanto que a secção mais baixa 
será 10 mm, 

A escolha da secção correspondendo ao seg- 
mento que, no diagrama, limita superiormente o 
domínio de cada tensão satisfaz o critério que 
se enunciou. Contudo este critério merece mais 
alguma reflexão. 

O critério significa que se escolherá a secção 
que permite produtos IP iguais aos que são 
suportados pelos condutores de 10 mm” quando 
se usa a tensão imediatamente superior, supondo 
constantes todos os outros parâmetros. 

Ora a expressão (7) é da forma 


Ka 
r + Kb 


1P = U*c . 


sms GD 


Se forem Uíc e Ux as duas tensões contíguas 
no diagrama, rz2 a resistência quilométrica do 
condutor de 10 mm? e r; a resistência quilomé- 
trica da secção si (mm?) que se pretende deter- 
minar, pode-se escrever, atendendo ao produto 
IP se manter constante, 


U?, 


U*ac 
ri + Kb 


=— É na s ga vã » IA 
r2 + Kb Ee) 


em que, para condutores em cobre, 


Kp= 0,3 ohm. quilómetro”! 


r2 == 1,72 ohm. quilómetro”! 
17,2 «14 
rr =“ ohm. quilómetro”! 
st 


Substituindo estes valores na expressão (11) 
obtém-se : 


(12) 


(1) Artigo 79.º do Projecto de Regulamento de Segu- 
rança de Linhas Eléctricas de Alta Tensão, Nov., 1962. 
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Para as tensões normalizadas, ou de uso mais 
frequente em Portugal, o Quadro 1 resume os 
resultados da aplicação da expressão (12). 


QUADRO 1 
Ulic U2Ze si 
(kV) (kV) (mm) 
[o 10 40 
10 15 28 
15 30 83 
30 60 83 


Em vez das secções obtidas pelo cálculo acima, 
e que não existem no mercado, adoptaram-se 
para definir o limite superior do domínio de 
cada tensão as secções comerciais que constam 
do Quadro 2. O limite inferior foi sempre defi- 
nido para a secção de 10 mm. 


QUADRO 2 

| Uic 51 
(kV) | (mm) 

6 35 

10 35 
15 70 
30 70 
60 


70 | 


Nas linhas de média tensão já não é cómoda 
nem tem muito interesse, a montagem das sec- 
ções de 50 a 70 mm?, As secções de 10 a 35 mm” 
são as de uso mais frequente. 


4, Há uma limitação para as potências trans- 
portáveis que importa não esquecer: a que re- 
sulta do aquecimento dos condutores. 

O parágrafo único do artigo 76.º do «Projecto 
de Regulamento de Segurança das Linhas Eléc- 
tricas de Alta Tensão», publicado em Novembro 
de 1962, recomenda que as intensidades da cor- 
rente em condutores de cobre não ultrapassem 
os valores indicados na segunda coluna do Qua- 
dro 3. A essas intensidades correspodem para 
cada tensão Uc, as potências que se indicam no 
mesmo quadro. 
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QUADRO 3 


Secção | Intensidade máxima 
admissível em regime 

| permanente 

(valor eficaz) 

(mm) (A) 6 
10 90 750 
35 175 1450 
70 280 2300 


Conclui-se que cada tensão será vantajosa- 
mente utilizada no domínio definido pela expres- 
são (7) mas a potência transportada nunca deve 
exceder os limites indicados no Quadro 3. Esta 
limitação justifica os segmentos paralelos ao eixo 
das abcissas que aparecem na Figura 1. 


5. O diagrama serve para definir zonas de 
utilização das tensões normalizadas, podendo, 
evidentemente haver interesse em examinar con- 
dições de transporte fora dos domínios que se 
consideraram na Figura 1. 

Para extrapolar os valores indicados na figura 
basta usar a expressão (7) 


10m Ufc 
r+x tgó 


IP = 


que permite concluir que o produto IP (ou qual- 
quer dos seus factores) cresce proporcionalmente 
com a queda de tensão m, e com o quadrado da 
tensão Uc. 

O aumento da secção reduz r e permite au- 
mentar 1IP mas em razão menor do que a pro- 
porcional, por se manter praticamente constante 
a reactância x. 

Embora se trate de relações elementares bem 
conhecidas, é vantajoso vê-las presentes ao uti- 
lizar o diagrama de Figura 1. 


6. A resistência quilométrica dos condutores 
que se tomou para estabelecer a Figura 1 está 
indicada no Quadro 4. 

O diagrama poderá ser aplicado a linhas cons- 
truídas em ligas de cobre ou alumínio, em alu- 
mínio-aço ou quaisquer outros materiais desde 
que a resistência r dos respectivos condutores 
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Potência activa limite P 
(em kW, para cos ( = 0,8) 


Tensão (kV) 
10 15 E 
1240 | 1850 | 3750 | 7 400 
2400 3650 7 300 14 600 
3900 5750 11 500 23 000 


seja análoga à dos condutores que figuram no 
diagrama. 


QUADRO 4 


Resistência unitá- 
ria (r) 


Secção do condutor em 
cobre duro (s) 


| 
(mm?) (2. km—!) 
O 172 
35 0,493 
70 | 0,246 


7. A tensão de 60 kV não é considerada uma 
média tensão. Não pode esquecer-se que supor 
nula a admitância transversal das linhas a 60 kV 
introduz erros apreciáveis no cálculo da queda 
de tensão. Contudo estendeu-se o diagrama da 
Figura 1 até aos 60 kV por ser essa a tensão 
normalizada imediatamente superior aos 30 kV 
ficando assim delineados os domínios de utiliza- 
ção das duas tensões contíguas. 

A figura refere, também, a tensão de 10 kV. 
Esta tensão, embora não mencionada no artigo 83.º 
do Decreto-Lei n.º 43 335, atrás referido, é usada 
com frequência e está normalizada pela C. E. I. 
(C. E. 1.-38-1954). 


8. Outros critérios são conhecidos e devem 
ser aplicados, paralelamente ao que se expõe, 
para definir o domínio de tensões que é vanta- 
joso a um determinado transporte de energia, 
Cita-se, como exemplo, o valor dado pela fórmula 
de Still corrigida, que é sempre uma indicação 
com interesse para os transportes em média 
tensão. 

O diagrama que se representa na Figura 1 
tem-se revelado, porém, um instrumento de tra- 


balho útil. 
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eia e TT 


JAYME DA COSTA, L” 


MECÂNICA E ELECTRICIDADE EM TODAS AS APLICAÇÕES 


PORTO LISBOA LUANDA 
Praça da Batalha, 12 Rua dos Correeiros, 14 Cx. Postal 5174 
Telef. 22812 (PPC) Telef. 327035 (pPc) Telef. 4743 


É aço, “e 


à É aretd 


Material fabricado nas nossas oficinas, pronto para expedição 


Fabricação de aparelhagem eléctrica de alta e baixa tensão 
Equipamentos electromecáânicos completos para centrais e fábricas 
Postos de transformação tipo monobloco 

Quadros blindados, capsulados e em armário 

firmaduras de iluminação de todos os tipos 

Contadores eléctricos e instrumentos de medida 

Telemedida e telecomando 

Instalações elevatórias de água 

Instalações de saneamento e esgotos 

Ferramentas de corte 

Motores Diesel industriais e maritimos 

Fogões eléctricos 

Frigoríficos 

Ventilação e condicionamento de ar 

Aiparelhagem de Raios X e eleciromedicina 
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ESTAMPARIAS TINTURARIAS 
LAVANDARIAS 


PARA OS VOSSOS TECIDOS TÉCNICOS: 
TRANSPORTADORES, MANCHONS, LAPPING, WOLFRIES, 
FLANELAS, CANEVAS, ETC,, 


EM FIBRAS NATURAIS OU SINTÉTICAS 
CONSULTEM: 


FANAFEL 


FABRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 


APART.: 9 
OVAR ESCRITÓRIO E FÁBRICA 


TELEFS : 52093-52094 
TELEG.: FELTROS PORTUGAL ESTRADA DE S. JOÃO 
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TECNICA Nº 359 


INTERNATIONAL SOCIETY OF ROCK MECHANICS 
1 st INTERNATIONAL CONGRESS 


LISBON / PORTUGAL — 25 TH SEPTEMBER 1 ST OCTOBER 1966 


The purpose of the Society's First Internatio- 
nal Congress is to provide an opportunity for 
the presentation and comparison of the ideas 
and experience on rock mechanics and to stimu- 
late their discussion by specialists from all the 
fields of science and technology interested in 
this new branch of mechanics. 


The Congress is intended to cover all the 
fundamental aspects of Rock Mechanics with 
interest for civil engineering, mining enginee- 
ring and petroleum engineering. It has been 
recognized that no case histories should be pre- 
sented, 

The following eight themes were considered : 


1. Exploration of rock masses: geologic explora- 
tion; boring and other exploration procedures; 
geophysical exploration method. 

2. Description of rock masses with a view to their 
physical and mechanical behaviour: texture; struc- 
ture; fabrics; anisotropy; heterogeneity ; joints; 
faults ; statistical methods ; quality indices; clas- 
sification of rocks and rock masses. 

3. Properties of rocks and rocks masses: physical 
properties; deformation and rupture under static 
and dynamic, tectonic or artificial actions; in- 
fluence of water ; alteration. 

4. Residual stresses in rock masses: measuring 
procedures ; stress distributions, their relations 
with tectonics, 

5. Comminution: physical bases of comminu- 
tion; breaking; drilling; blasting; crushing; 
grinding; abrasion. 

6. Natural and excavated slopes: stress distribu- 
tions; stability; stabilization of slopes; observa- 
tion of the behaviour. 

7. Underground excavations and deep boring : stress 
and strain distributions ; stability of escavations 
and boreholes ; linings; strengthening; support; 
subsidence; observation of the behaviour. 

8. Behaviour of rock masses as structural founda- 
tions: bearing capacity of rock masses; settle- 
ments ; observation of the behaviour. 


All correspondence concerning the Congress 
should be adressed to: 


O objectivo do 1.º Congresso Internacional da 
Sociedade é dar uma oportunidade de apresentar 
e comparar as ideias e a experiência adquiridas 
no que respeita a Mecânica das Rochas e esti- 
mular a sua discussão por especialistas dos dife- 
rentes domínios da ciência e da tecnologia que 
estejam interessados neste novo ramo da mecá- 
nica. 

Pretende-se neste 1.º Congresso cobrir total- 
mente os aspectos fundamentais da Mecânica das 
Rochas com interesse para a engenharia civil, 
de minas e do petróleo. Decidiu-se que não de- 
vem ser apresentadas descrições de casos parti- 
culares. 

Foram considerados os oito temas seguintes 


1. Prospecção de maciços rochosos: reconhecimento 
geológico ; sondagens e outras operações de pros- 
pecção; métodos geofísicos de prospecção. 

2. Descrição das rochas e dos maciços rochosos do 
ponto de vista do seu comportamento físico e mecânico ; 
textura; estrutura; «fabrics»; anisotropia; hete- 
rogeneidade; diáclases; falhas; métodos estatis- 
ticos; índices de qualidade; classificação das ro- 
chas e dos maciços rochosos. 

3. Propriedades das rochas e dos maciços rochosos : 
propriedade físicas; deformação e ruptura sob 
acções estáticas e dinâmicas, tectónicas ou artifi- 
ciais; influência da água; alteração. 

4. Tensões residuais nos maciços rochosos: processos 
de medida ; estados de tensão, suas relações com 
a tectónica. 

5. Desagregação de rochas: bases físicas da desa- 
gregação; fragmentação; perfuração; explosão; 
pulverização ; abrasão. 

6. Taludes naturais e de excavação : estados de ten- 
são; estabilidade; estabilização dos taludes ; 
observação do comportamento. 

7, Escavações subterrâneas e sondagens profundas: 
estados de tensão e deformação ; estabilidade das 
escavações e dos furos de sondagem ; revesti- 
mento; consolidação ; suporte; aluímento; obser- 
vação do comportamento. 

8. Comportamento dos maciços rochosos como funda- 
ções ; capacidade de suporte dos maciços ; assenta- 
mento; observação do comportamento. 


Toda a correspondência respeitante ao Con- 
gresso deve ser dirigida a: 


Secretariado do 1.º Congresso Internacional de Mecânica das Rochas 
Laboratório Nacional de Engenharia Civil 


Lisboa 5 — Portugal 


Telegraphic adress: LNEC Lisboa. 
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TÉCNICA Nº 355 


C. D. U. 536.252 


ESTUDO EXPERIMENTAL DA TRANSMISSÃO DE CALOR 
POR CONVECÇÃO FORÇADA NUMA CONDUTA DE 
SECÇÃO ANULAR 


por A. DE OLIVEIRA FALCÃO 


RESUMO 


Apresentam-se os resultados do estudo experimental da 
transmissão de calor por convecção forçada de água numa 
conduta de secção anular, com Ds/D, = 2,07 e aquecimento 
na parede interior, Para 9620 < Re <, 36100, a equação 
Nu = 0,0265 ReS Pril3 representa os resultados obtidos 
com um erro inferior a 109. 


1. INTRODUÇÃO 


Neste artigo apresentam-se os resultados do 
estudo experimental da transmissão de calor por 
convecção forçada de água, numa conduta de 
secção anular, com libertação uniforme de calor 
na parede interior. 

O imperfeito conhecimento do mecanismo da 
turbulência e as enormes dificuldades analíticas 
que comporta o estudo da convecção em regime 
turbulento restringem considerâvelmente o inte- 
resse prático dos resultados obtidos por via teó- 
rica, mesmo para casos relativamente simples. 

Tem-se recorrido, por isso, ao estudo experi- 
mental, para um grande número de configura- 
ções geométricas e condições aos limites, espe- 
cialmente os de maior interesse nas aplicações em 
engenharia, e procurou-se ampliar o campo de 
aplicação dos resultados por intermédio da aná- 
lise dimensional. 

É considerável a quantidade de resultados pu- 
blicados sobre a convecção forçada em regime 
turbulento em condutas de secção anular, em- 
bora seja pouco satisfatória a concordância duma 
grande parte deles. Efectivamente, este tipo de 
geometria tem larga aplicação no campo da en- 
genharia, nomeadamente em permutadores de 
calor e reactores nucleares. 
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Assistente do 1.5. T. 
Núcleo de Estudos de Engenharia Mecânica do LI, A.C, 


SINOPSYS 


The results of an experimental work on heattransfer to 
water fowing in an annulus are reported, for Ds| Dr = 
2,07 and heat transfer across the inner wall, For 9620 <, 
lte <, 96100, the results are correlated by theequation Nu == 
0,0265 ReS Pril3, with a deviation less than 109. 


2. APLICAÇÃO DA ANÁLISE DIMENSIONAL 


A análise dimensional é um método que se 
aplica, duma maneira geral, aos fenómenos que 
possam ser descritos por equações dimensional- 
mente homogéneas e baseia-se em que uma re- 
lação funcional, dimensionalmente homogénea, 
entre m quantidades expressas num sistema de n 
unidades fundamentais pode ser reduzida a uma 
relação entre (m —n) produtos adimensionais. 
Esta proposição constitui o chamado teorema 
dos = ou de Buckingham [1]. 
A grande importância deste teorema consiste 
em permitir reduzir considerâvelmente o número 
de variáveis do problema e ampliar a informa- 
ção obtida a partir dum número limitado de re- 
sultados experimentais. 
No nosso caso particular, consideramos que a 
quantidade de calor q transmitida por unidade 
de tempo e por unidade de área da parede 
[QL-* T-!] é função das seguintes oito va- 
riáveis: | 
Ds: — D; = De — diferença entre os diâmetros da 
secção anular [L] 

D:/D; — relação entre os diâmetros 

v — velocidade média do escoamento [LT”!] 

At — diferença entre a temperatura da parede 
e a temperatura média do fluido [O] 
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k — condutibilidade térmica [QL-! T-160-1] 

cp — calor específico a pressão constante 
[QM-! 0-1] 

o — massa específica [ML 


. 


u — viscosidade [ML-! T”!] 


Ao escolhermos um sistema de cinco unida- 
des fundamentais em que, além da unidade de 
temperatura O), se inclui a unidade de quanti- 
dade de calor Q, implicitamente admitimos que 
é desprezável a quantidade de energia mecá- 
nica que se transforma em energia calorífica e 
vice-versa. Mais particularmente, supomos que 
a energia dissipada em atritos é muito pequena 
comparada com a quantidade de calor permutado, 

Em segundo lugar, consideramos que o com- 
primento da conduta é suficientemente grande 
para que os efeitos das extremidades se não fa- 
çam sentir, ou seja, admitimos que os perfis da 
velocidade e da temperatura estão completamente 
estabelecidos. 

Aplicando os métodos da análise dimensional 
|1], podemos reduzir a relação entre as nove 
grandezas referidas a uma relação entre quatro 
grupos adimensionais : 


Nu = f (Re, Pr, D;/D,) (1) 
sendo 
Nu = q De Es h De o número de Nusselt 
St k k 
Re= "É ia o número de Reynolds 
a 
H Cp , 
Pr== —— o número de Prandtl 


As propriedades físicas do fluido k, cp, pe u 
são funções da temperatura e, portanto, variam 
de ponto para ponto. Esta variação é particular- 
mente acentuada para a viscosidade, no caso de 
muitos líquidos, quando as variações de tempe- 
ratura são consideráveis. Para ter em conta esta 
circunstância, é usual fazer figurar duas variáveis 
representando a viscosidade do fluido, referidas 
a duas temperaturas definidas. Neste caso, pas- 
samos a ter uma relação entre cinco grupos adi- 
mensionais : 


Nu = f (Re, Pr, DyDi,4/p'). (2) 
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Note-se que o comprimento De = D: — Di, 
que figura em Nu e Re, é o diâmetro equivalente, 
tal como é usualmente definido em hidrodinã- 
mica: De ==4A/P, sendo A a área da secção e P 
o perímetro molhado. Estudos experimentais 
mostraram que as leis das perdas de carga para 
condutas de secção circular de diâmetro D se apli- 
cam aproximadamente a outras formas de secção 
desde que se substitua D por De. Efectivamente, 
desde que seja constante a relação A/P, é a 
mesma a relação entre o gradiente da pressão na 
direcção axial e a força tangencial média por 
unidade de área da parede. Comparando com a 
secção circular, virá D==De. O conceito de 
diâmetro equivalente foi posteriormente aplicado 
ao estudo da transmissão de calor por convecção 
forçada, com resultados relativamente satisfató- 
rios. Alguns autores sugerem a utilização dum 
«diâmetro térmico equivalente», definido como 
sendo Dim = 4 A/Pimn, sendo Pm a parte do 
perímetro da secção ao longo do qual se realiza 
a transmissão de calor. 


3. CORRELAÇÕES EXPERIMENTAIS 
ANTERIORES 


Um número considerável de investigadores 
procurou estabelecer, com base em dados expe- 
rimentais, a forma da função (1) ou da fun- 
ção (2), tendo sido realizados ensaios com vários 
fluidos, mais frequentemente o ar ou a água. 
Apresentamos seguidamente algumas das corre- 
lações mais conhecidas, para o caso de trans- 
missão de calor na parede interior (!) (exclui-se 
o caso de fluidos com Pr muito baixo, como é o 
caso dos metais líquidos) : 


Monrad e Pelton [3] 
Nu = 0,020 Re"* Pr!” (Ds/Dy)"* (3) 


Davis [4] 
Nu = 0,038 Reº* Prló (u/uç)!4 
+ 0,2 4 
E o 1) (D:/D1)º!º (4) 
D; ; 


McMillen e Larson [5] 
Nu = 0,023 Re?º Prl!/º (D;/D,)"*º (5) 


(1) Um estudo mais completo pode ser encontrado 
num trabalho de Pickering (2). 
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Walger [6] 
Nu = 0,021 Re? Pri (us/pw)t (Da/Di)** (6) 


Os índices w e f indicam que as propriedades 
físicas são consideradas à temperatura tw da 
parede e à temperatura tr = (tw + tm)/2. Às 
restantes propriedades são tomadas à tempera- 
tura média tm, definida como sendo 


tm=([tepvds)/( [co ev ds) 


Note-se que os mesmos expoentes 0,8 e 1/3 
de Re e Pr figuram em todas as expressões apre- 
sentadas e que coincidem com os das expressões 
clássicas de Colburn e de Sieder e Tate, obtidas 
para tubos de secção circular [7]. 


4. DESCRIÇÃO DA INSTALAÇÃO 


Na nossa instalação utilizámos água destilada 
escoando-se numa conduta de secção anular 
(fig. 1), com D:= 20,3 mm e D;=9,8 mm, limi- 
tada exteriormente por um tubo de vidro e in- 
teriormente por um tubo de aço Cr-Ni, com a 
espessura de 0,14 mm e comprimento L = 500 
mm. O diâmetro equivalente é De = D; — D; = 
10,5 mm e D:/D; = 2,07. 

O tubo de aço Cr-Ni tem uma resistência de 
aproximadamente 0,10 ohm e é percorrido por 
uma corrente eléctrica fornecida por um trans- 
formador, de intensidade variável até 250 Amp, 
a que corresponde uma potência máxima dis- 
sipada de cerca de 6 kw. A determinação da 
potência fez-se medindo a intensidade da cor- 
rente e a tensão aos terminais do tubo, com o 
auxílio dum multímetro de 1,5 "/o de precisão. 

Para a determinação da temperatura da parede 
do tubo, foi utilizado um termopar de Ferro-Cons- 
tantan com fios de 1 mm” de secção, que ficou 
encostado à superfície interior do tubo, podendo 
ter deslocamentos angulares e longitudinais. A 
temperatura do termopar foi medida por um gal- 
vanómetro com uma precisão de 1ºC, estando a 
soldadura fria mantida a temperatura constante 
num banho de gelo fundente. 

Para se obter a temperatura tw da superfície 
exterior do tubo, foi necessário entrar com a 
diferença entre as superfícies interior e exterior 
do mesmo, dada por q a*/2 k”, sendo q a quanti- 
dade de calor libertado por unidade de volume 
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do tubo e por unidade de tempo, a a sua espes- 
sura e k a condutibilidade térmica do material 
do tubo. Esta diferença tem o valor máximo de 
1,2ºC para os ensaios efectuados. 

As temperaturas do fluido à entrada e à saída 
da conduta anular foram medidas por termóme- 
tros de mercúrio, com divisões de 0,2ºC. 

A circulação de água é assegurada por uma 
bomba de carretos, accionada por um motor 
eléctrico de velocidade variável. O líquido, à 
saída da conduta de ensaios, é conduzido a uma 
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caixa de descarga sobreelevada donde, por gra- 
vidade, desce directamente para o reservatório 
ou atravessa primeiramente um balão de vidro, 
conforme a posição duma válvula colocada junto 
ao reservatório. O caudal é calculado medindo 
o tempo necessário para encher o balão, antes e 
depois de cada ensaio, tomando-se a média dos 
dois valores. 

No interior do reservatório existe uma serpen- 
tina para refrigeração, destinada a manter a esta- 


caixa de 
descarga 


Ci lreservatório 
Ce em 
ES serpentina de 


= ! refrigeração 


bilidade térmica do sistema para os vários valo- 
res da temperatura da água e da potência liber- 
tada. 


5. RESULTADOS DOS ENSAIOS 


Os ensaios foram realizados para uma secção 
situada a x==350 mm da entrada da conduta 
anular. A relação x/De==33 é suficientemente 
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grande para que os efeitos de entrada se não fa- 
çam sentir, pelo menos para o regime turbulento. 
A temperatura média tm da água na referida 
secção foi calculada a partir das temperaturas 
de entrada te e de saída ts: 


tm = te + (x/L) (ts — te) 


O coeficiente de transmissão de calor h foi 
calculado com base na diferença de temperatu- 
ras (tw Cie tm). 


As propriedades físicas da água foram toma- 
das para a temperatura tm [8]. 

Foram realizados 31 ensaios, para Re variável 
de 1230 a 36100, apresentando-se, na tabela 
junta, os respectivos resultados. 

Na fig. 3 estão representados, gráficamente, 
os valores de Nu/Pr !" em função de Re. 

Para Re > 9,62 >< 10º, os pontos obtidos afas- 
tam-se menos que 10"/y da recta: 
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Nu = 0,0265 Re"$ Prlf 
e, para Re > 15,3>< 10º, o desvio em relação à 
mesma recta é inferior a 3º/,. 


Para o valor utilizado D>/D; = 2,07, os coefi- 
ciente de proporcionalidade das relação (3) a (6) 


em 13º/, ao valor 0,023 usualmente tomado para 
tubos de secção circular [7]. 

Para 2,91>< 10º < Re < 8,53 >< 10º, os pontos 
obtidos ficam abaixo da recta referida. 

O regime laminar está representado por 4 
pontos, com 1,23> 10º < Re < 1,81 >< 10º, ali- 
nhados aproximadamente segundo uma recta com 


tomam os valores 0,0294 (Monrad e Pelton), 
0,0430 (Davis), 0,0320 (McMillen e Larson) e 
0,0292 (Walger). O coeficiente 0,0265 obtido é 
inferior a estes em 10"/) a 62!/, mas é superior 
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o coeficiente angular de 0,35, valor que não se 
afasta muito do expoente de Re encontrado 
por vários autores para tubos e condutas anula- 
res [9], 


Q G tm 
kw g.s-1 OC 
3,94 366 66,4 
3,94 345 | 66,7 
3,78 333 67,9 
3,48 267 67,3 
2,97 226 67,1 
3,92 253 49,6 
2,47 163 66,6 
3,92 246 49,5 
3,73 224 49,3 
3,17 201 49,1 
2,84 201 | 41,6 
2,82 176 48,3 
2,35 155 47,5 
2,83 156 42,5 
2,19 131 47,1 
MIS 131 40,7 
1,78 107 46,4 

2,36 101 40,5 1 
1,40 98,7 36,3 
1,05 78,0 35,5 
1,17 79,3 33,6 
0,955 73,4 33,0 
0,965 63,7 33,2 
0,774 64,2 31,8 
0,689 52,1 32,9 
0,608 41,3 31,1 
0,455 35,0 30,6 
0,447 33,3 31,1 
0,447 28,7 31,0 
0,448 25,8 | 31,5 
0,448 22,4 31,9 
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Pr Re »<10-3 Nu 
2,70 36,1 161,0 
2,69 34,2 159,1 
| 2,64 33,7 | 151,3 
| 267 26,7 130,1 
2,68 22,5 110,6 
3,56 19,4 111,4 
2,69 19,1 988 
3,56 18,8 105,4 
3,57 17,1 96,5 
3,59 15,3 87,7 
| 4,22 13,4 90,8 
3,63 13,3 74,2 
3,70 11,5 65,1 
4,12 10,5 67,9 
3,73 9,62 58,7 
| 4,28 8,53 52,4 
3,77 7,78 46,4 
4,28 6,54 40,9 
4,68 5,96 40,5 
4,76 4,65 33,6 
4,99 4,54 31,9 
5,05 4,15 29,6 
5,03 3,62 26,8 
5,21 3,52 27,6 
5,11 2,91 23,7 
5,29 2,23 22,7 
5,36 | 1,87 19,9 
5,27 | 1,81 16,0 
5,30 1,55 15,3 
5,24 1,41 14,6 
5,20 | 1,23 13,8 
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DIMENSIONAMENTO DE FUNDAÇÕES DE MÁQUINAS 


II 


FUNDAÇÕES DE MÁQUINAS ALTERNATIVAS 


RESUMO 


Este artigo é o segundo de uma série subordinada ao 
titulo geral «Dimensionamento de fundação de máquinas», 

Apresentam-se as condições gerais a que devem obedecer 
as fundações de máquinas alternativas; estudam-se as goli- 
citações dinâmicas impostas pelo funcionamento destas má- 
quinas às fundações. Com base no artigo anterior da 
presente série, apresenta-se como aplicação numérica o 
cálculo da fundação de um motor Diesel de 2 cilindros: 
Enumeram-se os métodos para limitar a amplitude de 
vibração das fundações, limitam-se os respectivos campos 
de acção, indicam-se normas de cálculo para amortecedores 
de molas, borracha ou cortiça e finalmente apresenta-se o 
cúleulo de uma fundação com amortecedores mecânicos 
de molas, 


1 — INDICAÇÕES GERAIS 


1.1 — Amplitude máxima de vibração das 
fundações 


Na generalidade das máquinas alternativas, 
caracterizadas pela transformação de movimen- 
tos de vai-vem em movimentos de rotação, de- 
senvolvem-se esforços não equilibrados que se 
transmitem às fundações. Dado que a frequência 
destas máquinas é em geral baixa, as vibrações 
podem ter influência considerável na estabilidade 
de estruturas ou de fundações adjacentes, pelo 
que se torna de primeira necessidade uma aná- 
lise, tanto quanto possível exacta, da amplitude 
das vibrações desenvolvidas. 

É difícil estabelecer um limite para a ampli- 
tude das vibrações. Há casos em que esta atinge 
valores de 0,4 a 0,5 mm sem inconvenientes sen- 
síveis. Noutros casos, porém, vibrações de ampli- 
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SINOPSYS 


This is the second article of a series entitled aDesign 
of engine foundations», 

The general conditions to which the foundations of 
reciprocating engines shall obey are here presented. Dy. 
namic strains imposed to the foundation by the running 
are studied. It is presented, as a numeric application, the 
calculation of a 2-cylinder Diesel motor, which is based 
upon the former article of this series. Methods used to limit 
the amplitude of vibration of the foundation are here quoted 
their fields of action being limited, Some rules are indica- 
ted for the caleulation of springs, rubber or cork absorbers 
and the compulation of a foundation with mechanica 
spring load-absorbers is presented, 


tudes inferiores, com frequências baixas, trans- 
mitem-se através do terreno, danificando estru- 
turas colocadas a dezenas de metros. 

Tomaremos como valor limite da amplitude 
nas fundações de máquinas alternativas 


A = 0,20 a 0,30 mm 


1.2 — Especificações necessárias ao estudo 
da fundação 


Os elementos necessários, a fornecer pelo fa- 
bricante, são os seguintes: 


1 — Velocidade e potência da máquina. 

2 — Características, grandeza e ponto de aplica- 
ção das solicitações dinâmicas devidas ao 
funcionamento da máquina, ou elementos 
necessários para o respectivo cálculo. 

3 — Distribuição das cargas estáticas impostas 
pela máquina. 
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4 — Dimensão e forma da base da máquina. 

5 — Localização de aberturas na fundação des- 
tinadas a parafusos de amarração, tuba- 
gem, etc. 


1.3 — Material de fundação 


Em máquinas Diesel a carga total, estática e 
dinâmica, transmitida à fundação, é tal que a 
tensão na face superior da fundação não exceda 


3 a 5 kg/cm? 


Em máquinas de cilindros horizontais este va- 
lor é ainda menor. 

Em qualquer dos casos, as tensões no bloco 
de fundação são pequenas, limitando-se a sua 
análise a pontos singulares, correspondentes a 
furos ou irregularidades de forma. 

Emprega-se correntemente, em fundações de 
máquinas alternativas, betão pobre de resistência 
à compressão aos 28 dias 


Res = 100 kg'emº 


1.4 — Normas para o dimensionamento das 
fundações 


Geralmente as fundações deste tipo de máqui- 
nas são constituídas por blocos maciços, com as 
aberturas e furação necessárias. Devido a essa 
forma, é lícito aplicar-se no seu estudo a teoria 
de vibrações dum sólido assente sobre um meio 
elástico, referida em [1]. 

A condição básica para o dimensionamento do 
bloco de fundação é a de que as suas dimensões 
sejam as mínimas para que a amplitude das vi- 
brações não exceda o limite admissível. 

Para que o assentamento da fundação seja 
uniforme, o centro de gravidade do conjunto 
máquina-bloco e o centro da área da base da 
fundação devem estar na mesma vertical. Caso 
isto não seja possível, a excentricidade não deve 
em caso nenhum exceder 3/o da dimensão me- 
nor da base da fundação, pelas razões apontadas 
em [1]. 

O bloco de fundação deve ficar sempre com- 
pletamente separado das estruturas adjacentes 
(sapatas, paredes, pavimentos, etc.). Não deve 
servir de fundação a outros elementos ou má- 
quinas; no caso de não ser possível evitar a 
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colocação de elementos estruturais secundários 
sobre o bloco de fundação da máquina, é indis- 
pensável atenuar as vibrações daqueles por inter- 
posição de placas de cortiça, borracha, feltro ou 
outros materiais isoladores. 

No caso de se instalarem várias máquinas na 
mesma oficina a distâncias semelhantes às di- 
mensões dos blocos de fundação, é aconselhável, 
em solos brandos, colocar a fundação de máqui- 
nas semelhantes sobre um estrado comum, de 
rigidez tal que as suas deformações sejam infe- 
riores às amplitudes de vibração dos blocos. 
Como o cálculo das vibrações de um conjunto 
destes é muito complicado, é corrente considerar 
o estrado dividido em blocos associados à fun- 
dação de cada máquina e calcular as vibrações 
de cada fundação em separado, considerando 
uma amplitude limite superior ao normal em 
cerca de 20 a 25º/,. 

As várias chumaceiras de um veio, bem como 
as chumaceiras de máquinas directamente aco- 
pladas a máquinas alternativas devem ser mon- 
tadas sobre o mesmo bloco de fundação, para 
evitar deformações no veio. 

Quanto maior a área duma fundação, mais 
elevadas serão as respectivas frequências pró- 
prias. Este facto é muito importante, pois per- 
mite, no caso das máquinas em estudo, que as 
frequências próprias das fundações sejam muito 
superiores às frequências forçadas pela máquina. 
A massa duma fundação deve pois distribuir-se 
de modo a que seja mínimo o momento de inér- 
cia em relação ao eixo horizontal que passa pelo 
centro da base. Portanto, deve limitar-se ao mi- 
nimo a altura dos blocos de fundação. 


2 — SOLICITAÇÕES DINÂMICAS IMPOSTAS 
PELAS MÁQUINAS ALTERNATIVAS 


As vibrações forçadas em fundações de má- 
quinas alternativas devem-se principalmente às 
forças de inércia, não equilibradas, dos elemen- 
tos móveis: êmbolo, haste, biela e manivela. 


2.1 — Máquina de um cilindro 


Consideremos o esquema de máquina alterna- 
tiva representada na fig. 1. A posição do sis- 
tema em qualquer instante é definida pelo ângulo 
de rotação da manivela, 2, 
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Fig. 1 — Esquema da máquina alternativa simples 


A — êmbolo 
B — haste 

C — biela 

D — manivela 


Qualquer elemento i do sistema dá origem a 
forças de inércia Pi de componentes segundo os 


eixos xx e zz respectivamente iguais a 
Pxi = m;.xXi 
tt 
Pi = mi . Zi 
sendo mi a massa do elemento i 
As forças de inércia totais serão 


di 
“ 


i 


Sendo w a velocidade angular, suposta cons- 
tante, da manivela, será 


FP, d” Zi 
X=w—— 
d 9? 
ti dº Z 
Zi= w' 
d 9º 
e vem 
d” Zi 
Pu wu 2 ml — 
i d 7 2 
+ ” d” Zi 
Po == w' 2 mi —— 
i do 2 


Sendo xt, X2, X3, Z1, Z3, Z3, as coordenadas 
do centro de gravidade respectivamente da mani- 
vela, conjunto êmbolo-haste e biela, e Mt, Mi, 
Ms, as respectivas massas, será = 
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dº xa dº xs d” xs 
mid 
Pe = wº (M: FE + Ms do? + Ms do? ) 
d? x4 d? zo d' 73 
sam 
P,=w* (M; da + Ma do + Ma FT ) 


É possível concentrar a massa do sistema ape- 
nas em 2 pontos, a e b, indicados na fig. 2, o 
que simplifica as expressões de Px e P, 


s| 


Fig 2 - Localização dos pontos a e b onde se considera 
concentrada toda a massa do sistema 


Demonstra-se ser: [3] 


Pe =R wº. Ma sen wt 
P, == Rwº [ (Ma + Mb) cos wt + 
(2) 


cos wt] +... 


as 


+ «Mp (1 + =! cos 2 wt — MA 


sendo as massas equivalente Ma e Mp dadas por: 


Ri Ls 
Ma ==— M — M 
E 
Mp = Ma +! Mh 
R 
q == — 
L 


designando por R, Rá, L, Ly, Ls, as grandezas 
indicadas na fig. 2. 

É de notar que as forças de inércia segundo 
o eixo Z Z contêm uma harmónica fundamental, 
que corresponde à força de inércia primária e 
harmónicas superiores, que correspondem às 
forças de inércia secundárias. Nas eg. (2) despre- 
zam-se os termos de ordem superior, por o seu 
coeficiente ser rapidamente decrescente. 

Instalando contrapesos de massa M'a sobre o 
veio, é possível anular as forças devidas a Ms. 
Assim, se M'« for colocado de forma a que o 
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seu raio vector forme um ângulo 7 com o raio 
vector de Ma (ou seja, no prolongamento da 
manivela) e a uma distância Ide O (fig. 2), para 
que as forças de inérçia dos elementos que rodam 
se anulem, (P, == 0) deverá verificar-se: 


ou seja (3) 


Para anular as forças de inércia primárias, na 
direcção da embol haste (zz), deverá ser: 


Ma + Mb — Ma =0 


ou (4) 
É + Ma - Mi — M' =0 
R E” 


Se os valores | e M, obedecerem só a uma 
das condições, na expressão de P, não se anulam 
os termos secundários, contendo cos2wt e cos4wt, 
e o valor de Px será: 

Pe == Rwº. Mp sen wt 

Normalmente Mp —> Ma e portanto a condi- 
ção (4) leva a um aumento das forças de inércia 
segundo o eixo XX, 


2.2 — Máquinas de vários cilindros 


Seja uma máquina com 1 cilindros em associa- 
ção linear sendo Pk a desfasagem angular entre 
a manivela do cilindro K e a do 1.º cilindro, de 
referência. Consideremos um sistema de eixos 
segundo indica a fig. 3. 

Desprezando as forças de inércia secundárias 
de coeficiente de ordem igual ou superior a 2º, 
será : 


e 
M cilindro 


Peilindro z 


Fig. 3 — Esquema de associação linear de cilindros 
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n 
Pe = wº 2 Re Mak sen (wt + fx) 
ad | 


P; = wº Fa Rk [ (Mar + Mp) cos (wt + Px) + 
+ Mp «k cos 2 (wt + 4] (5) 


Estas forças dão lugar a momentos segundo 
os 3 eixos, de valor 


n n 
Mi = » Pak o lyk ; My = 2 Pak o lx; 
k=1 k=1 
n 
Mz sua ça Pxk . lyk (6) 
Se os cilindros forem idênticos, será: 
E 
Pa = Rwº Ma , 2 sen (wt + fx) 
o ; 
P, = [(M, + Mb) R as (wt + Po) + 
| à 
+ Mp x A cos 2 tw +) 
] 


Portanto, para equilibrar as forças de inércia 
primárias, deverá ser: 
n 


Z cos (wt + fi) = 
k=1 


Para anular as harmónicas secundárias deverá 
ser: 


n 
O; > sen (wt + Ex) =0 


Z cos2(wt+fj)=0 
k= 
Para equilibrar os momentos devidos às forças 
primárias, será : 
lx cos (wt + Ex) = O 
lx sen (wt -H Cx) = 0 


lyx sen (wt + 2x) — O 


mis ris ria 


Vejamos alguns casos correntes, considerando 
cilindros iguais e desprezando as harmónicas su- 
periores. 


2.2.1, — Máquina de 2 cilindros verticais, montados asst- 
métricamente sobre a fundação 


1.º Caso: manivelas em fase: 


Será fi=0 fe=27 
Pa = Pa = Rwº. Ma . sen wt 
P,1 == P,) = Rw? (Ma + Mb) cos wt 
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Fig. 4 — Esquema de máquina com 2 cilindros verticais 


Py; = 2 RwºM, sen wt 

P, = 2 Rwº (Ma + Mb) cos wt 
Ms = 2 P,1 (1 + 21,) 

My = 2 Pw lz 

M; =P, 1+21,) 


Este tipo de máquina é fortemente desequili- 
brado, dando origem a vibrações intensas. 


2.º caso: manivelas desfasadas de 90º 


Será agora: 


Pa == Rwº M, sen wt 


Px == Rwº Ma sen (w + E = Rwº Ma cos wt 


P, = Rwº (Ma + Mb) cos wt 
P,2 == (Me + Mb) cos (mt + dl = — 


= — Rwº (Ma + Mb) sen wt 
ou seja 


v==V2 Rwº Ma sen (w + =) 
P, =V2 Rwº (Ma + Mb) cos (m + 2) 


As componentes resultantes são V2 vezes 
maiores que as de cada cilindro. 
Os momentos serão neste caso: 


Mes as Pi (1 + ly) + Ps) ly 
My = (Px + Pxs) 1; 
M, = Pu (1 + 15) + Pwzly 


Este caso é o mais corrente em máquinas de 
2 cilindros. 


3.º caso: manivelas desfasadas de 180º 


será 
Pa = 


Ps; = 0 


Bs == O 


P, == 0 
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e ainda 


Ma == Pa l: My =0; M, = Pu l 


2.2.2 Máquinas de 2 cilindros horizontais 


É caso corrente em compressores, com as ma- 
nivelas desfasadas de 90º. As forças de inércia 
são dadas pelas expressões indicadas para o 
2.º caso, trocando x por z nas equações. 

É ainda: 

Mex a Ps; lyi o P;o ly2 
Ms = P; lx + Pk | 
Mz = Pa ly1 — Py2 lya 


2.2.3 — Máquinas de 3 cilindros verticais 


Normalmente é 


Bb =0; fPa= 120º; [bs = 240º 
Vem: 
Px =0; P, = 0 
My = P,1 (21 + ly) + Po (1 + 1,) + Pz3 ly 
My = 0 


M, = Ps (21 + ly) + Ps? (1 + ly) + Pws ly 


2.2.4 — Máquinas de 4 cilindros verticais 


E: fi=0; Ba = 180º; Es = 180º; Bs = 360º 
Anulam-se tanto as forças de inércia como os 


momentos. 


2.2.5 — Máquinas de 6 cilindros verticais 


Neste caso, é geralmente : 


2 


ti=0 i=— fB=b= e Bs = 27% 


As forças de inércia primárias são nulas. Os 
momentos são: 


M«a=V3 Pr 1I M=0 M=V3 Pl 


Os seus valores são pequenos pelo que em 
fundações de máquinas de 8 cilindros não é em 
regra necessário o cálculo de vibrações forçadas. 


3 — SOLICITAÇÕES IMPOSTAS POR TRANS- 
MISSÃO POR CORREIAS 


É muito corrente o caso de máquinas alterna- 
tivas associadas por meio de correias a outros 
tipos de máquinas. 

O movimento da correia dá lugar a esforços 
sobre as chumaceiras das máquinas e conse- 
quentemente sobre as fundações. 

Seja a ligação por correia representada em 
esquema na fig. 5. 


rodo mbtora 


Fig. 5 — Esquema de transmissão por correias 


Sendo T; e T; os esforços sobre as correias, a 
força P, exercida sobre cada chumaceira será 


P: = Ti + Ts 
A força periférica transmitida pela correia é 
Pp = Ti— Ts 


Sendo W a potência da máquina e » a velo- 
cidade da correia, é 


' W 
Pp = — 
v 
como 
2 TE 
v = Rw =— NR 
60 
em que: 


N — velocidade de rotação da máquina, r.p.m. 
R — raio da roda não motora 


vem : 


Pr = = 9,55 (8) 


7 — menor ângulo de contacto entre correias 
na roda motora 
wu — coeficiente de atrito entre correia e roda 


obtém-se : 

W 1 

Tj== 9,55 ——* (9) 
N R 3 e!” 
W pu 

De E os 
N R T e” 
W 1I+He” 

Fp=s 9,55 — “——— (10) 
NR 7 e! 


A direcção da força Pt é a da linha que une 
os eixos das rodas. 

Se esta direcção fizer um ângulo y com a 
horizontal (xx), a componente horizontal de P; 
será: 

Pix = Pt. cos 7 


A componente vertical de P, é em geral des- 
prezável, por ser muito pequena em comparação 
com o peso da máquina. 

Se as duas rodas estiverem montadas sobre a 
mesma fundação, as forças desenvolvidas pela 
transmissão são forças internas, sem influência 
sobre o equilíbrio da fundação. 


4. APLICAÇÃO NUMÉRICA: CÁLCULO DA 
FUNDAÇÃO DE UM MOTOR DIESEL 


De acordo com os elementos de Mecânica dos 
Solos e Teoria das Vibrações expostos no artigo 
anterior [1], segue-se o cálculo do bloco de fun- 
dação de uma máquina alternativa muito cor- 
rente, o motor Diesel de cilindros verticais. 


a) Características do motor ; 


Potência: W = 27 Cv 

Cilindros: 2, verticais 

Velocidade de regime: 1500 r.p.m. 
Área da base: 0,70>< 0,75 m 
Peso: 600 kg. 


b) Características do terreno : 


É constituído por argilas brandas, sendo apro- 
ximadamente : 


— tensão de segurança à compressão, em re- 
gime estático: 1,5 kg/cm” 
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— coeficiente de compressão elástica uniforme : 
Cu = 3,0 kg/cm” 

— coeficiente de compressão elástica não uni- 
forme: C; == 6,0 kg/cm” 

— coeficiente de corte elástico uniforme : C, = 
= 1,5 kg/cmº. 


c) Solicitações transmitidas pelo motor 


Não sendo conhecidas com exactidão as carac- 
terísticas necessárias à determinação de P, e P, 
(2.1), admitimos ser: 


Mi = M:=M) 
R = 10 em 


É fácil relacionar, dentro de certa aproxima- 
ção, a potência do motor com as massas equiva- 
lentes (2.1). O movimento do êmbolo e da haste 
será da forma: 

x = À sen wt 


O trabalho desenvolvido por este movimento 
em cada ciclo tem o valor: 


— Ms Aº wº sen wt 
O 


cos wt dt=— 2 M, Rº wº 


referindo-se os simbolos à notaçãs empregada 
até aqui. 
A potência é então: 


W = 50M:;Rº Wº = 27 Cv 


As massas equivalentes são neste caso 
Ma=M =Mi =mM 
Mb un 1,5 Ma 

A força desenvolvida por cada cilindro, de 


acordo com as equações (2) tem componentes 
horizontal e vertical iguais respectivamente a 


P'yx == 205 sen wt (kg) 
P', = 512 cos wt (kg) 


Admitindo um desfasamento entre manivelas 
de 90º, as solicitações resultantes são: 


Py = 2 Py = 290 sen wt (kg) 
P, = V2 P, = 720 sen wt (kg) 
My == 256 sen wt (kgm) 
My, = 205 sen wt (kgm) 
M == 102 sen wt (kgm) 
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d) Forma da fundação 


Dada a conveniência de que as frequências 
próprias da fundação sejam sensivelmente supe- 
riores à velocidade de regime da máquina, que 
é elevada, adopta-se um bloco de pequena altura 
de forma indicada na fig. 6. 


Fig 6 — Fundação de motor Diesel 


O peso da fundação será: 


P=1.815 kg 


Cota do centro de gravidade do conjunto má- 
quina-fundação : 


e) Amplitudes de oscilação 


Consideremos, em 1.2 análise, os diversos 
modos de vibração como independentes. Será 
então : 

Frequência própria das oscilações verticais : 

pt — CUA  3>x10x20x<1,7>9,8 

m 2415 
== 41.400; f, = 203 rads” ! = 32,5 Hz 


Amplitude das oscilações verticais : 
Pz 720 9,8 


À ui Da E agi 
“ m(n?-f2) 2415(41.400 — 24.800) 


= 1,75 10"* m == 0,175 mm 


valor inferior ao limite admissível. 

Oscilações angulares transversais 

Sendo o momento dinâmico transversal 
My == 256 sennt (kgm) 


2RN 
60 


tw 


= 50 1 


A . [4 * - 
e a frequência própria de rotação 


: CGI-PZ Ci 


f.;* visto ser PZ « Cyl 
Wo Wo 
em que: 
= 6 kg/cmº 
00 = 170º : 
[= ec = 8.25< 10" em" 
12 
2415 = 30º " 
Wo =——— —— = 2.210 kgems” 
981 
Vem: 


fo? = 223.000; fy = 470 rads”! = 75 Hz 


a amplitude de oscilação angular é 
eua “M o 256 o 
Wo(m?—f;?) 22,1 223.000 — 24.800) 
= 5,78>< 107º rad 


As 


a amplitude de oscilação vertical correspondente 
nos bordos da fundação é: 


A = 5,78>< 107º x 0,85 m == 0,049 mm 


No caso mais desfavorável de sobreposição 
de A, e Azl, a amplitude vertical resultante 


será 
A= A, + Àz1= 0,225 mm 


Valor que está dentro do limite admissível 
(0,20 a 0,30 mm). 


Oscilações horizontais, segundo xx ou yy 


Sendo a respectiva frequência própria: 


CA 15x 10º x 3,40 9,81 
A Ei DE 4 fai 
m 2.415 


= 20.800 ; 
ft = 144 rads”* = 23 Hz 


de valor muito próximo da velocidade de regime 
do motor (25 Hz), torna-se necessário aumentar 
a rigidez de corte da fundação, o que se pode 
conseguir embebendo a fundação no terreno e 
compactando-o fortemente na zona próxima do 
bloco. Procedendo assim, o coeficiente de corte 
elástico uniforme Cz aumenta considerávelmente, 
bem como a frequência própria horizontal. 
Admitindo-se um acréscimo de Cxz para 2,5 
vezes o seu valor inicial, o que está de acordo 
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“A . A 
com a experiência [2], o novo valor da frequên- 
cia própria será : 


E, = 52.000; f,= 228 rads”! = 36,3 Hz 


que já é sensivelmente superior à velocidade de 
regime. 

A amplitude das oscilações horizontais, será: 
DO . 290 x 9,81 o 
“m(Ê—w?) 2415 (52.000 — 24.800) — 

— 4,35 x 10º m = 0,044 mm 


Ax 


f) Tensões sobre o terreno 


Tensão estática vertical: 


2415 ' 
fg = — = 0,071 kg/cm” 
3,4 >< 10º 


Pressão dinâmica vertical: 
sa = ACu = 0,067 kg/cm* 
Pressão vertical total: 


= + ca = 0,138 kg/cm? 


Verifica-se assim que embora a tensão sobre 
o terreno tenha um valor muito reduzido, a 
amplitude das vibrações verticais é da ordem de 
grandeza do limite admissível, o que indica cla- 
ramente que em problemas desta natureza Oo 
objectivo fundamental do dimensionamento das 
fundações é limitar a amplitude das oscilações 
transmitidas ao terreno. 


5 — MÉTODOS PARA LIMITAR A AMPLI- 
TUDE DE VIBRAÇÕES DAS FUNDAÇÕES 


5.1 — Equilíbrio das forças de inércia, por 
meio de contrapesos 


Como se viu em 2.1, é possível, por meio de 
contrapesos : 


— equilibrar por completo a força de inércia 
normal à direcção da haste do pistão e parte 
da componente na direcção da haste; 

— ou equilibrar por completo a harmónica fun- 
damental de componente segundo a direcção 
da haste do pistão, do que resulta um acrés- 
cimo para a componente normal àquela: 
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Normalmente emprega-se o 1.º método, por 
serem menores tanto os contrapesos como as 
tensões na máquina. 

Depois de instalada a máquina, se se verificar 
que se desenvolvem vibrações exageradas, recor- 
re-se ao método que as limita, dependente da 
natureza das vibrações. 

O equilíbrio por meio de contrapesos causa 
apenas uma curta interrupção no funcionamento 
da máquina. 


5.2 — Estabilização química dos solos 


Aplica-se especialmente em solos arenosos. 
A estabilização aumenta a rigidez da base e 
portanto as frequências próprias de vibração, o 
que é vantajoso, como se viu, em máquinas deste 
tipo. Pode ser aplicada rápidamente, sem grande 
interrupção no funcionamento da máquina. 

Se as vibrações forem sobretudo angulares, 
pode limitar-se a estabilização aos bordos da 
fundação. 

Em regra basta levar a estabilização à profun- 
didade de 1 a 2 metros. 


5.3 — Processos estruturais 


São os mais dispendiosos e requerem uma 
interrupção de funcionamento das máquinas ge- 
ralmente longa. Empregam-se portanto apenas 
quando os métodos anteriores não dêem resultado. 

Consistem: na cravação de estacas embebidas 
na fundação, que lhe aumentam considerável- 
mente a rigidez; na ligação da fundação a placas 
especiais, destinadas a aumentar a rigidez do 
conjunto; [2] no emprego de massas ligadas à 
fundação, que funcionam como amortizadores 
[1] [6] [7] e se revelam muito eficientes nos 
casos em que não há grandes variações na ve- 
locidade de regime da máquina. 


5.4 — Amortecedores mecânicos de molas 


Empregam-se amortecedores de molas sempre 
que é necessário limitar as amplitudes de vibra- 
ção muito aquém dos valores normais. 

As molas são de fácil instalação, eficientes e 
pouco dispendiosas. São de muitos tipos, prote- 
gidas geralmente por cortiça, borracha ou feltro. 

As características das molas a empregar (diã- 
metro e número de espiras) dependem dos re- 
sultados de cálculo dinâmico. 
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Conforme o tipo da máquina e a sua veloci- 
dade, existem 2 dispositivos de montagem das 
molas: de suporte ou de suspensão. 

O 1.º caso aplica-se em fundações de máqui- 
nas de alta velocidade (superior a 300 r.p.m.) 
bem equilibradas, em que não é necessário um 
bloco de fundação de grandes dimensões. As 
molas são então instaladas directamente sobre a 
fundação. Sobre elas assenta um quadro metá- 
lico rígido, sobre o qual está montada a máquina, 
conforme indica em esquema a fig. 7. 


bloco de 


fundocõo 


Fig. 7 — Esquema de fundação de molas de suporte 


Pode ser necessário aumentar a massa sobre 
as molas. Nesse caso, ligado ao quadro metálico 
dispõe-se um bloco de fundação ligeiro, assente 
sobre as molas, colocadas por sua vez sobre o 
bloco de fundação principal. 

Para a instalação dum sistema deste tipo, a 
ordem de operações é a seguinte: Betona-se o 
bloco inferior da fundação, com uma altura ge- 
ralmente de 1 a 2 metros, dependente do tipo 
de máquina e das propriedades do solo; depois 
do betão endurecer, reveste-se a face superior 
com borracha sintética, aglomerado betuminoso 
ou qualquer material análogo; colocam-se os su- 
portes das molas; assenta-se o bloco superior ; 
montam-se as molas e regulam-se; finalmente, 
monta-se a máquina. 

As molas de suspensão empregam-se em má- 
quinas de baixa velocidade em que se torna 
necessária uma fundação pesada sobre as molas, 
a fim de evitar a colocação destas a um nível 
inferior ao do pavimento, com os inconvenientes 
de difícil regulação e substituição. 

Neste caso os parafusos que ligam as molas 
ao apoio da máquina são de grandes dimensões. 
O conjunto tem o aspecto indicado na fig. 8. 

Como as molas são especialmente empregadas 
em máquina de cilindros verticais, a análise di- 
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NANA 


bloco superior montado 


NINA 


obre o quadro metalico 


ANNAN 77 


£< 


8 
N 
A 


ro m 
s1e 


—- e — msm em — cem 


mola 


porafuso de fixação 


NAN Ti 


ou de suspensão 


N 


À 
DAN, 


cortão asfaltico 


ou borracha 


Fig. 8 — Esquema de fundação de molas de suspensão 


nâmica reduz-se ao estudo de vibrações verticais 
de um sistema com 2 graus de liberdade. 

Seja então um sistema de massas concentradas 
mi e ms, ligadas por molas de rigidez we c, 
actuado por uma solicitação periódica P sen wt, 
conforme se indica na fig. 9. 


Psenwtl 


C; 


C4 


Fig. 9 — Modelo de 2 graus de liberdade, representativo 
do caso em estudo 


As equações que regem o movimento do sis- 
tema são, desprezando o amortecimento: 
tt 
m z:+czi—c(z2—zi)=0 
(11) 
m:z2+ ca (22—2z1) =P sen wt 


sendo: 
ms — massa dos elementos colocados acima 
das molas (máquina e fundação) 
my — massa da fundação colocada abaixo das 


molas 
zt,  z3—elongação vertical das massas mem» 
ci =cu À — coeficiente de rigidez do solo 
cu — coeficiente de compressão elástica uni- 


forme do solo 
A — área da fundação 
cz — coeficiente de rigidez do conjunto de 
molas, dado por: 
nd 


— 6 LIZ 
8n Dº ) 


Cs = 
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m —n.º de molas 

n —n.º de espiras em cada mola 

D — diâmetro das espiras 

d — diâmetro da mola 

G — Módulo de elasticidade transversal do 
material das molas 


Considerando apenas vibrações forçadas, será: 


zi= Assenot; Zzs= As sen ot 


| == o 
my À (1?) 


(1 +)? + ve —n? 


“ 
+ 


Às = :P (13) 
ms A (192) 
sendo: 
Ap Êo Ms o ca NR 
ma mi + ms mi 


A (o) =0t—(1+0) (BHE)RA (+ PE? 


Vejamos a relação entre as amplitudes A, e A» 
e a rigidez Cs» das molas. O valor da força P é 
proporcional a w”, ou seja 


P= tw? 


sendo y dependente das características da má- 
quina. 

Substituindo esta expressão em (13) e divi- 
dindo por ww”, obtém-se: 
ig é E e OM 
ml 1-(1+v) E +): E?) 


sendo : 
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Na ausência de molas, a amplitude tem o 
valor : 
- 1 
ne A. : 


2. (15) 
m (1+v) (:72—1) 


Portanto o coeficiente de absorção das molas 
será : 
A 1-0+) Er? 
Ai (1+pv) &*—1) 


7 — 


Vejamos então a relação entre £ e n: 


Para f<<w,ou seja sei — 0 


vem: n— o9 

Quer dizer: se a frequência própria dos ele- 
mentos colocados acima das molas, supondo 
estas assentes sobre um meio infinitamente ri- 
gido, for muito inferior à frequência forçada, 
a eficiência da mola é muito grande. 

Suponhamos agora que É -—> co, o que cor- 
responde a um grande valor de f neste caso 
n— 1, ou seja, a mola não tem qualquer efi- 
ciência, 

Apresenta-se na fig. 10 a variação de n com É, 
na qual se torna patente que as molas são efi- 
cientes para [1 | > 1, o que sucede para uma 
faixa de valores É compreendida entre O e És 


tal que 
RE VA o] (1 +- H) DP =] (17) 
2(1 +04) (621) 


Fig. 10 — Relação entre x e 


Parai >%,|7n|]<1,ou seja, o emprego de 
molas é inconveniente, pois aumenta a ampli- 
tude das vibrações. 


Portanto, para que as molas sejam eficientes, 
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é necessário que a frequência própria da massa 

montada sobre elas seja sensivelmente inferior 

à frequência forçada. Esta frequência própria 

depende da rigidez das molas e da massa ma. 
Se se fixar o valor de n, obtém-se: 


1+)(n—1) (67? —1) 
5.5 — Tapetes de borracha e cortiça 


São de emprego muito corrente, como siste- 
mas de limitação da amplitude das vibrações. 

O seu comportamento é inteiramente seme- 
lhante ao das molas, pelo que o respectivo cál- 
culo é regido pelas normas apresentadas em 5.4. 

No caso de apoios de borracha, o cálculo 
apresenta certas dificuldades pela acentuada não- 
-linearidade deste material, ou seja, da variação 
dos módulos de elasticidade com as condições 
de carga e forma de apoio. 

É de notar ainda que a borracha necessita ser 
protegida da acção de agentes corrosivos e tende 
a alterar-se com o tempo, o que levou à sua 
substituição por diversos produtos sintéticos. 

Indicam-se nas referências [7] e [10] da bi- 
bliografia, os elementos necessários ao cálculo 
estático e dinâmico de apoios de borracha. 

As propriedades mecânicas dos aglomerados 
de cortiça empregados correntemente no isola- 
mento de vibrações são descritas na publica- 
ção [11] da lista bibliográfica. 


5.6 — Aplicação numérica: cálculo da funda- 
ção de um compressor vertical 


a) Dados: Trata-se de um compressor vertical 
de 120 kw. As vibrações devem ser reduzidas 
ao mínimo, pois existem próximo aparelhos 
de precisão. O solo é caracterizado por um 
coeficiente de compressão elástica uniforme 
Cu = 2>< 10º ton/m, 

A amplitude das vibrações não deve exceder 
0,03 mm. 

A velocidade do compressor é de 480 r.p.m. 
A força de inércia vertical é P, = 2,6 ton. 


b) Amplitude das vibrações 


Suponhamos que a relação entre a frequência 
própria da fundação e a frequência forçada é 2, 
Para que fosse: 

Az = 0,03 x 10º m 
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O peso da fundação deveria ser: 


2,6 < 9,8 
its cana, — = 115 ton 
0,03>x 107º X 3x 2,5 >< 10º 
e a área da fundação seria: 
2W 4x 2,5x 10ºx 115 . 
fa E = 57,5 m 


Cu 2 x 10º 9,81 


Estes valores são excessivos, dada a pequena 
potência da máquina, pelo que é conveniente 
recorrer à utilização de molas. 

Escolheu-se a fundação indicada em esquema 


na fig. 8, sendo: 


— área da fundação: A = 12,5 m” 

— peso da fundação abaixo das molas (bloco 
inferior): 21,0 ton 

— peso do bloco superior e máquina: 35,0 ton 

— coeficiente de rigidez da base: c = cu A = 
= 2 x 10º x 12,5 = 25,0 >< 10º ton/m 


Vem portanto : 


21,0 

m=—— = 2,15 ton seg? m 
9,81 
35,0 

m: = ——— = 3,56 ton seg” m”! 
9,81 


— frequência própria do sistema, suposto sem 


molas : 
c ; 
= =4,38>x 10 seg”? 
ms + ma 
* f 2 
= = 1,76 
ul 
= E qa 1,65 
mi 


Na ausência de molas, a amplitude vertical 


seria: 
2,6 
As = = em 
(2,12 1 3,56) (4,38 - 2,5) 10º 


= 0,25 > 10"“m == 0,25 mm 


O coeficiente de absorção das molas deverá 
então ser: 
As (0,25 
Pav 0,03 


f = 8,4 
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Suponhamos que 4 =— 10 . Obtém-se, da 
eq.. (18): 
CC 1T—(1+1,65)>x1,75, 
(1 + 255) (= 70=3) (1,75 —1) 


.a 
 “ 


= (1,165 


O valor da frequência própria da massa as- 
sente sobre as molas será então 


f;* = E? w? = 414 seg 
e a rigidez total das molas: 
ca = f,* ms = 1.480 ton'm 


A rigidez de cada mola é então 


m 

mi —n.º de molas 

seja ny — 16; vem: 
— 1480 


c=>— = 920 ton/m 
16 


Por outro lado, é: 
. 4+G 
8D'n 
G — módulo de elasticidade transversal =— 


7,5 > 10º ton'm? 
Sejam molas de 5 espiras. Vem: 


7,5><10º da 1 A 
92,0 = — — , 1 sexos É 
D 5 D' 
d* 
= =4,9 x 107 
Dº 
Seja : 
d=2,5><10""m; 
vem: 


D=9,3 cm 
c) Análise da tensão nas molas 


A capacidade de carga de uma mola é dada 
por 
ed? £ 
Po SE; a = (19) 
8 D 
t — Tensão de segurança à torsão: 
40 = 10º ton/m? 


Obtém-se para este caso: 


Po = 2,64 ton 
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Para determinar a carga que actua de facto 
sobre cada mola é necessário calcular a ampli- 
tude das vibrações forçadas. Será: 


A (nº) = 6,25>< 10º — (1+1,65) 
(0,414 10º + 4,48>< 10º) x 2,5 >< 10º + 
+ (1 + 1,65)><0,41 >< 10º 4,38 x< 10º = 
== 21,0 x 10º 


vem então, da eg. (13): 


(1 +1,65)><4380 +1,65><414 - 2500 
3,56 x 21 >< 10º 
>< 2,64 = 0,34>< 10“ m 


As >< 


A carga dinâmica sobre as molas é então: 
Pa = 0,34>x< 10º X 1,48>< 10º = 0,5 ton 


A carga total será 


Pe == id al À = 2,22 ton < 2,64 ton (limite 
10 admissível) 
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